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ABSTRACT 


An analysis is made of termperature distributions in a slab of finite 
thickness and infinite extent, one surface of which is perfectly insulated 
and the other surface of which is exposed to a fluid the temperature of 
which varies in a specified way. Two cases are considered. In the first 
case, the temperature of the fluid is suddenly changed and maintained at 
the new value. In the second case the temperature of the fluid increases 
linearly with time. From the solutions of these problems, curve sheets 
are developed which permit evaluating each of eight different quantities 
of engineering significance. These curve sh2zets cover more cases of 
interest and are easier to use than similar curves which have appeared in 


the report literature. 
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I. INTRODUCTION 


The problem of convective heat transfer by a fluid flowing through a 
thin-wall pipe with a well insulated exterior surface is of considerable 
importance, part:iicularly in connection with assessing the effects of rapid 
operating transients in nuclear power systems. Several years ago contractors 
engaged in this sort of work began using a set of several sheets of curves 
which facilitate these calculations. The writer has been unable to learn 
the original source of these curves or even to find them anywhere in the 
open literature. However, they are to be found in such Wocwiente oe Re fer- 
ences 2 and 3. 

The theory upon which these curves are based is well known and will be 
developed, for reference purposes, in Section II of this thesis. It re- 
places the thin-wall tube by a flat plate which extends infinitely in the 
y and z directions and is of finite thickness a in the x direction. At 
the surface x = 0 it is presumed that there is a perfect insulating layer 
and at the surface x = a the plate is exposed to a fluid whose temperature 
varies with time in a specified manner; the coefficient of surface heat 
transfer between the fluid and the plate is a specified constant h. The 
initial temperature distribution in the plate is given as a function of x. 
The problem is to determine the subsequent temperature distribution in the 
plate as a function of x and the time t. It is further presumed that the 
appropriate physical properties of the plate are constants. 

This problem can be reduced to that of solving a certain linear partial 
differential equation subject to certain initial and boundary conditions. 


To simplify and systematize the analysis it is convenient to introduce and 





deal with several dimensionless ratios. There is more than one form in 
which the solut:ion may be represented but we will obtain and develop a 
particular and convenient form. 

In the case of a pipe the wall erictrecese? which is an appreciable 
fraction of its radius, cylindrical geometry should be employed and an addi- 
tional geometric parameter, ratio of wall thickness to radius, must be intro- 
duced. The solution in this case is more complicated than that which will. 
be developed here and which forms the basis for the useful curves referred 
to above. However, in most pec eiC a cases the wall thickness is indeed 
sufficiently small compared to radius to justify using the flat-plate geo- 
metry which wi]l be employed here. No attempt is made in this thesis to 
explore the differences between results using flat-plate geometry and re- 
sults using cylindrical geometry. 

Two particular cases arise frequently in practice. In each case it is 
assumed that the fluid and the material of the pipe are at the same tempera- 
ture prior to the time t = 0. In one case, referred to as "step", it is 
assumed that the fluid suddenly assumes and maintains a new and different 
temperature starting at time t = 0. In the other case, referred to as 
"ramp'' or "linear", it is assumed that starting at time t = 0 the temperature 
of the fluid Reneee linearly with time. By superposition of results, the 
cases of piecewise linear changes in fluid temperature including sudden 
"jumps" in temperature can be accommodated, if convenient solutions are 
available for the two cases described. 

The useful curves referred to above give significant information 
relating to these two cases. However, those sets of curves which have 
come to hand are considerably less useful than would be desirable because 


they do not provide a sufficiently broad range of significant parameter 





values to represent physical situations actually encountered and they require 
some rather drastic interpolation and extrapolation. Furthermore, certain 
integral properties of the solutions are useful for certain types of stress 
calculation and the curve sets previously available have not provided these 
properties. It should also be noted that some of the information contained 
in the curves which are developed could also be obtained from the well known 


' specifically from Figures 7 and 10 of Reference 4. 


edeusler charts,' 

Accordingly, it has been a major purpose of this thesis to reconstruct 
the curve sets, including a greater variation in significant parameters and 
to augment the collection of curves by adding several new curves likely to 
prove useful in engineering practice. The curves developed were checked and 
verified to be correct by comparison to previously published curves. The 
only exceptions to this were the AT, and AT. cases for which there were 
no comparison curves or data available. 

The development, including the differential equation and boundary and 
initial conditions, leading to the production of the desired curves is 
presented in Section II.A. of this thesis. In Section II.B. the theory of 
the step fluid temperature change is developed and that of the linear fluid 
temperature change is developed in Section II.C. A discussion of the compu- 
tational method used to obtain the curve sets is presented in Section II.D. 
Additionally, the computer program has many comments inserted to aid in 
understanding it. A brief discussion of how to use the curve sets is pre- 
sented in the first part of Appendix A, following which the curves them- 
selves are presented. Further instruction is given in References 2 and 3. 


The person who is engaged in heat transfer calculations involving slabs or 


thin wall pipes can easily develop further methods as required. 





Inasmuch as it is likely that in nuclear power plant calculations there 
may be need to perform calculations relating to convective heat transfer 
of the type described above on a massive scale, involving hundreds of differ- 
ent parameter sets and hundreds of transient specifications, it is of in- 
terest to consider ways in which human intervention is not required to "read" 
curves. Of course, a computer may be used to calculate each of these possi- 
bly hundreds of thousands of cases using the exact theory developed herein. 
However, since the calculations are, as will be seen, quite demanding, there 
is reason to seek an approximation which will have sufficient accuracy for 
most purposes but which will provide significant speed-up of the calcula-- 
tions. For this purpose, in Section III hereof, a seek to determine the 
coefficients of a quickly-evaluated approximation in such a manner that the 
error will be limited to small values over the entire significant ranges 
of all parameters. 

Figure 1 illustrates the problem which will be dealt with. 


/ | | 
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We will find a convenient dimensionless temperature T(x,t)/T, for the 
two previously mentioned types of fluid temperature change. The step 


temperature change, where T, = 0 for time t < O and T. = F = const. for 


t > 0, will be considered first in Section II.B. The linear temperature 


a O for t < O and Ty = Et for t > O will then be con- 


sidered in Section II.C. Here E = const. is in degrees per unit time. 


change, where T 


Actual solutions will be cast in the dimensionless form T(x,t)/T, which will 
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be an explicit function of x or x/a and Fourier number and an implicit 
function of Biot number. 

Each of the two general equations giving T = T(x,t) will then be employed 
to provide eight particular properties most likely to be of use to the 
designer. Curves representing these properties in a convenient way will be 
obtained and exhibited in Appendix A. These results include those which 
‘have appeared in the curve sets of References 2 and 3 but with a fuller 
coverage of significant parameter values. Additionally some new properties 
are presented that are intended to assist in performing evaluations according 
to the rules of the ANSI B31.7 Code for Nuclear Power Piping and a forth- 
coming edition of Section III, Nuclear Pressure Vessels, of the ASME Boiler 


and Pressure Vessel Code. 
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II. THEORY AND DEVELOPMENT OF EQUATIONS FOR HEAT TRANSFER 


A. GENERAL THEORY 


It is desired to find T(x,t) such that it satisfies the unsteady heat 


conduction equation in one dimension, namely 

eee 

ee OL 

apie eps (1) 
The initial and boundary conditions of the problem are: 


1) T(x,0) = 0 (slab initially at uniform temperature throughout, 
and we take this to be zero without loss of 


generality) 
2) ot = 0 (perfect insulation at x = 0) 
OX =O 
3) a = hf -T (a; | (continuity of heat transfer at 
Oe ol Kt fluid-solid interface) 


B. STEP FLUID TEMPERATURE CHANGE 


Here TT. = F = const. for t > 0. The solution to this problem is 


Pound to be 


ed 
1%) - ae, Bese ” (2) 


where: 
a) P,, = M,~ G) =\, x 
b) An= My Va/A 
yen= Ma (A 





oennender: Paul J., Conduction Heat Transfer, p. 250, Addison Wesley, 


Jie) sye)e 
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d) C =4 sin M /(2M + sin 2M ) 
n n n n 
e) 2: denotes summation over n = 1,2,3,...00 


f) The numbers Mare an infinite number of distinct roots as 


seen in Figure 2 below, of the equation 


M tan M = N, (3) 


where Np is the Biot number 


N, = ha/k 
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Figure 2 


To verify that the governing differential equation is satisfied we find 









IT | a An 
ee Cr Po sin (By, X) © (4) 
wae ae Joti & 
PT 2 FS Cy Gf,” coslhnX) < 
me 
a- An i. 
* FE. Cp dn’ ces (b,x) © 


2 





The governing equation is then 


ey. = . Noe 
eS C, B,? cos( Gn x) eo” = FEL Cada’ cosCBnx) & 


But AB, = «Mis Lo? ty Therefore the differential equation is 
Satisfied. Next, we consider the first boundary condition which takes 
the form 

[= 3 Cy cosh, ¥) 
and we will show that this is true by actually determining the quantities 
Co These are similar to Fourier coefficients which must represent the 
constant value of unity. Even though the argument Gia x) is not simply 
a multiple of n, the customary technique of evaluation can still be used. 
Multiply through by cos( ee x)dx and integrate from 0 to a. On the left 


we get 


4 


. ie aS Coe x) dx = SG ie 


Din 


On the right we get 


aN 
ZiaGe .. ase, x) cos (fh, x) d x 


and it is necessary to distinguish between the cases where m =n andm#Zn. 


i E (Sin Ain) peer (Airmen) | 
oh errr re nein aE. 3 - 


[Aix = hp 14 int [Mr | 


E | ees Ie eos AN Kan TA i = Mr ALA bY Mj. 


For m gn, we find 


Sc 25 (fin %) GO) a= 


ip [ee 


which vanishes since MS and MA are distinct roots of M tan M = Na" 


For m = n, we find = 
A 7 Spee 
oe | ae _ ee 


e =a *CB, x)d 9 7 en ike cos” + il ai | Zz Mtr, J 


Thus, we have established that 
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os yy Sth 2. hii 
es s/n Mn ad Ce] | Zz Mn 


= / + . 
so that C_ + sin M /(2M_ sin 2M) 
as previously defined. 
That the second boundary condition is satisfied may be immediately 
verified by inspection of equation (4). Additionally, upon noting that 
fh 
y 1) Cer A) = Mn S17 Mp, = A Cos Mp 
X KK 
it is trivial to show that the third boundary condition is satisfied. 
Thus the governing equation and boundary conditions have been shown to be 


satisfied. Equation (2) is interesting and useful‘in its own right. 


However, we now proceed to derive formulas for eight particularly useful 
consequences of the relation. In these formulas the distance x no longer 
appears as a parameter and the dependence upon time is incorporated in the 


value of the Fourier number 
oa 
i Fee = 
(~) = Fourier number = % / 


which is simply related to the exponent in the general formula, viz: 
2s 
The eight cases we consider and the corresponding formulas are: 


1. Case 1. The average temperature of the slab is 





Cnn 
he ener 
Eve t SO THAD A, aL ][x-5) Se sine! | 
~ KR —@O ea CH pr, - 
sin Mn ef" C Mn @7 hn 
meet A-7, Cr sin tee fe >. ; s/f n 
a a Po é a ) 
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Pease une Neated suriace temperature is, by inspection, 
Ta, x - 
Ss MIs le Dy, Cp cos Mp er Fh 
3. Case 3. The insulated surface temperature is, by inspection, 


TO 2 os: Cr CTF 


rE 


4. Case 4. The difference between heated surface and average 


temperature is, by comparing Cases 1 and 2, 
pace, =) — 7. me ey 
e (4 


5. Case 5. The difference between average and insulated surface 


temperature is, by comparing Cases 1 and 3, | 
ee es cet (1 20M 
F ne My 


6, Case 6. The difference between heated and insulated surface 


temperature is, by comparing Cases 2 and 3, 
4 ~ oo ant ~~ ' 
ee CO) Ls BE Ph G eas M,) 


7. Case 7. The "Equivalent Linear Variation" is denoted by the 





symbol AT,- This quantity, which is defined by equation (5), is useful 
in determining what are known as secondary stresses due to radial thermal 


gradients in piping. 


~ 


A < (Xx - & ) BLS. dx 


ie 
7 aa 5 “pk Ta one fal x cos(fi, x )-8 = Salen 
tn 
2 7 [12 (1-cos Mn)-6M,, ip Mn | 


8. Case 8. The "Peak Component of Thermal Gradient" is denoted by 


the symbol AT,. This quantity, which is defined by equation (6), is 
useful in determining what are called peak stresses due to the nonlinearity 


of the radial thermal gradient in piping. 
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AT,/F = 


we { [7 (2, =)-Tay 2) £ fat], [ev )-T or -L [ar], ob /e - 


It should be remarked that the definitions of AT, and AT, given above 
are those which at the time of this writing have obtained the approval of 
the Working Group on Piping of the Design Su>committee on Section III of 
the ASME Boiler and Pressure Vessel Committes. It is probable that they 
will appear in a forthcoming edition of the Code, but several stages of 


approval remain before this can happen. 


C. LINEAR FLUID TEMPERATURE CHANGE 
Here T = Et for t > 0. The solution is aa by Duhamel's method 
which gives 
ne 
Memeo ny, <)(t-C *) 
gfx - & Sn eos CEnxi (ire (7) 


TOx,2*~) = xD, ® 


or 


(8) 


An 


roy x)= EAL a Cn 608 PaO" ~ 


since AL een ene 


To show that this solution satisfies the differential equation, we get 


from equation (8): 


An * 
OT. g r{S Bn Cn sin En (-€ vhe Je (oe (Gori Cin 510 By x)(I-e 0 *) 


aay 7 Mo 


ax in Nae 


aera notes of Professor John E. Brock 
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EOS : 
PT. ES tn cosh i-e 7) Leste, cos (Pax (I- onan *) 


Oxe An “SB 


and from _ ail 
ne 
oF = Eft Ser cos (@,x)C , | 


However, the constants C. were previously determined so as to satisfy the 
‘ —t 
equation fa. >: Gr eas Ce. x) o 


e-An'*) 
Thus we can write at ae Gaaes 5 CB, x) Cl 
and it is clear that we have satisfied the differential equation 
ieee 
ax * QA 


Now consider the boundary conditions. From equation (7) it is clear 


that T(x,0) = 0 and from our expression for at it is clear that | eC. 
IX : : IX X~O 


Finally, considering the third boundary condition, we have 


; ae *) 
oT | -~ E@ Ch Slo [yp Cie 
ae R= X Mp 

and finding T(a,t) from equation (7) gives ee 

| fin (ine *) 
AP Tard] = E4 of cy cos Minlle® 
— ei 

ae 
Paper Pe Ss 4 5 

- = is <1 Cp Sih) Ii, (i-€ ¢ iy ) ar I oo. 

H Fi, K Gn sie Min 


However, the last term in parenthesis is recognized as 
N/M tan M = 1 


so that the third boundary condition is also satisfied. The solution may 
now be expressed in the more desirable form 


eh 
TOM) 2 [- 5 Se cos C4, x) (I~ ¢ ) 
EX an - 





The eight special cases for the linear problem are, using the same 
notation as before: 


1. Case 1. 


IE 











0 

By) ot CTA) dys t]x-5 Sa om Ch, 0-2") 

ee 82 ez £, p 0 
1h 


= 1-5' C, sin M, (I-@*) 


M, FP. 
2. Case 2. By inspection 
ThE) pis Cy cos M, Ci-e) 
fe i p,, 


3. Case 3. By inspection 
ee 
TCC, *) [ass GEiCh= = ° 
EY p 
4, Case 4. From Cases 1 and 2 a 
BCGy Fy = (ay Cx) - §! SE anaes mw, (enh ic cosh \ 
Fe _ A 7 


5 
Bee Case 5. ‘From Cases 1 and 3 


Tay €#) ~ 16% *) 2 ss am ES eC. sin Bn 





[ere 
Gy Case 6, From Cases 2 and 3 


Peron) —Te*)_ = estce teensy 


Gees 


7. Case 7. Again using the eas eri nieton 


els 7 ie zo (x- =) rook Bl 


oe ry 


=>) Coin Or" I2 Cl- Cos M,)- 6M, Fl | 
(a | 


8. Case 8. Again 


EA 


D. COMPUTATIONAL SCHEME 


ee. vax [reo ey 448% |, @-TE, *)|- “2/41 all fe 


The curves, which are exhibited in Appendix A, were obtained by compu- 
tation on the IBM 360/67 computer at the Naval Postgraduate School. The 


programs which were used are presented in Appendix B. There should be no 
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difficulty in following these programs, which were written in FORTRAN IV 
language using double precision arithmetic; they are quite straightfor- 
ward implementations of the equations developed in the preceding sections 
of this thesis. | 

In covering the range peipacemeters invoi.ved, values of the Biot number 
were introduced in a dimensioned array while the Fourier numbers were 
generated by a DO-loop. A large number (41) of values of Fourier numbers 
were used; these are the abscissa values for the curves. A smaller number 
of values of Biot number meee Can ova: these correspond to the individual 
curves on each sheet. It was necessary to select pHese values by trial and 
error in each case to permit a reasonably attractive distribution of curves 
on each sheet. 


In calculating the roots of the equation M tan M = N., Newton's method 


B? 
was used; a root was declared to be established when the next approximation 
made less difference than ton The number of roots actually determined and 
the number of terms in the summations was thirty. No specific investigation 
was made of the effect of taking a larger or smaller number of terms; how- 
ever, some earlier work done by Professor Brock indicates that sufficient 
accuracy for graph plotting could probably have been obtained with only a 
few terms. 

A remark is called for with regard to the plotting routine. Although 
the calculations were performed with double precision arithmetic, the plot- 
ting routine, a local library subroutine called DRAW, accepts only single 
precision inputs. Therefore, as a final step prior to plotting,the appropri- 
ate arrays had to be changed from double to single precision format. 


As previously noted, the program should be easy to follow with the 


comments that have been inserted. If it is desired to change the number 


rg 





and/or value of the Biot numbers for any particular problem, the value of 
"RELL" for that problem and the particular BN(J) array would have to be 
changed. The width and/or interval spacing of the abscissae (log,)5@) 
can be changed by changing "KALL" (the number of points plotted for each 
Biot number) or "FANG" (the interval spacing for the abscissae). The starting 
point for this interval is set by the constant A in the card 

LOGTH(I) = A + (I-1)*FANG 
which appears in each particular problem. Many "WRITE" statements were 


included to provide a printout of data in addition to the graphs. 
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IIIT. DEVELOPMENT OF APPROXIMATE EQUATIONS 


A. THEORY 
As can be seen from the foregoing developments, the calculation of 
temperature at any given time in the interior of a slab which is convec- 
tively heated in the ways we have considered is a rather demanding task. 
Of course modern digital computers can perform such evaluations rapidly, 
but there are a large number of arithmetic operations involved in each such 
computation. Automated analysis ‘of complete piping systems subjected to a 
variety of thermal transient conditions might: require hundreds of thousands 
of such evaluations and in such a case it would be valuable to have a suffi- 
ciently accurate approximation which could be evaluated in a fraction of the 
time required for an evaluation using the theoretically exact equations. 
Thus we seek to find an accurate representation of T(x,t) in the form 
T = w(x,y,z) where x,y and z are conveniently chosen functions of x/a, Nee 
and @ respectively* and the function w is to be determined. We shall at- 
tempt to find the functions x, y and z so as to permit covering the values 


of x/a : d the many decades of N, and of 2, that we wish to cover and still 


B 


Peovide fairly uniform accuracy over the entire range. For the function w 


we can think of no more useful representation than a power series 


mM N P 
< y A in-1 net Pol 
1 >, > : map x y a eT) 
Inst fof Peal 


where the coefficients Aes are to be determined from a large number of 
point quadruples (WX 59%) obtained by use of the exact solution. A 


*Note that x which appears hereafter is a specified function of the para- 
meter (x/a). It should not be confused with the distance x which 


was used previously. 





Zk 





customary way of doing this is to determine the STG 50 as to minimize the 
sum of the squares of the errors (w, W(X, sy, 52,))5 where w is the approxi- 
mate solution and Ww. the known value from the exact solution. We now pro- 
ceed to establish the theory of determining the ernG to do this. Let the 
positive integers M, N, P , and Q be given where Q is the number of sets 


of four numbers x., Yeo Zz o and w, (i = 1, 2,...Q). Determine the A in 
i alt mnp 


equation (1) so that 


iMo- 


Cure wry a) 


is a minimum. The minimizing condition is 


- QM NP ; 
108 LSS ss oh al PL De boty ot RT 
COCSE zt frat asl f= = 


eee, Ms S = 1,2,...N; ¢ = 1,2,...P 


Et ets 2 | PA 7 5 aber sel | nl 
yy Pcs Vis aes 7 (2) 
Peel acl pst Ee { eet 


We now introduce the notation 


Ee | 
ea x ~{ ; 
aan = >, x7 yo a (Do not confuse this with 


i an earlier and different 


po use of the letter 8 ye 


ach gd =mt rel; K=n+s-l; L = p + t-l 


aa t-! s-l K-( 
and a ve Wh XK yo ee 
‘fae | 
Equation (2) becomes 
o N P 
Me Ag Pegs 7 Prs2 ™ 
, pie: pol bh) fp wick Poe 


ZZ 





Equation (3) is linear in the unknown coefficients A so that no 
mnp 


essential difficulty is involved in solving for them. If M, N, and P are 


small integers (say no larger than 4) a brute force type of arrangement 


and rearrangement of the equations will provide a form which can easily be 


solved by computer to yield the desired quantities. However, if the M, N, 


and P are larger integers, some sublety may be required in devising a 


computer progrem to perform the necessary manipulations. 


The problem of finding the correct values of the coefficients one be- 


comes one of: 


a) 


b) 


c) 


d) 


e) 


Selecting the three appropriate functions x, y, and 2 to’ be 
used, 


Selecting values of M, N, and P. Larger values should give 
move accurate results at the expense of more computing time 
and larger storage requirements. 

7 B 0 aput 
Generating values of Cone and ee the computer. 


Solving for the an from equation (3) by linear methods. 


Checking the approximate solution against the exact solution. 


As a simple example of this process consider the case where M = N = P = 2. 


This gives the following system. 


b, Bs Pes Boe can sae Boy Coe rw = M0 


one pe ee Bo i2. S13 Co22 br2? Arye Xie 


Bist Biz Baa br22 Car E232 A yet Vion 
B 133 F222 2.2.3 Porz P23 pie DG. 

Bort Cs) Bay co ee ce 

C 12 B 322. P5235 Ao Yee 


Symmetric 


Baz B 332 Aan Oo: Ze | 
6333 ee: ce Be 


<= 


Z 





The nD will provide an approximate solution over the range of data 


points used in the generation of pe and Oe ae: We used the functions 


x = cos(x/a), y = @ /(1+6)), and z = N,/ Ct +N.) and covered the ranges 


O<x/a £1, .001 < @ < 100., and .01 < N.< 1000000. 


B, COMPUTATIONAL SCHEME 
The first program,'Generation of Betas and Gammas," works as follows: 
1. The Biot numbers and Fourier numbers needed to generate the 
data points are established as for the heat transfer curves. 
2. 30 roots of M tan M = Na are calculated. 
3. The exact formula derived in Section II is used to seqeteee 
144 data points (i.e., Q = 144). 
4. The values of a and Nene are then calculated using their 
defining equations. | 
The beta and gamma matrices are now assembled from their known element:s 
and the SnD founds: by solving the system 
Le IA Le! 
for the unknown value of sae Then a program, hereinafter called "Check 


Program", uses the theory of Section III.A. to quickly generate the approxi- 


mate solutions. Values were printed for x/a QO and x/a = 1 and compared 


} 


to the results of the curves for x = 0 and x = a to verify the approxima- 
tion. This program displays the values of the ap for both the step and 
ramp cases and demonstrates how the approximate solution is used on the 
computer. Tne programs are shown in Appendix B and are self explanatory 
with the included comments. 


lone SMIS program of Prof. Giles Cantin was used to do this. 
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However, the checks were so unsatisfactory that it must be concluded 
that a programming error remains undetected or that our choice of the 


functions x, y, and z was very unwise. 
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IV. RECOMMENDATIONS FOR FURTHER WORK 


The approximate solutions generated by the "Check Program" in Appendix 

B were in such error that no specific recommendations can be made. It 
appears that either there was an error in the generation of the matrix 
elements and therefore in the evaluation of the ae or that the selection 
of the functions x, y, and z was unwise. Time limitations do not permit 
investigating this further. An attempt to correct the approximation should 
include; 

1. Review the "Generation of Betas and Gammas'' program to insure 
that it in fact does properly generate the values of ies and 2 awe 

2. Rewrite this same program so as to 


a) eliminate the argument INDEX from the arrays of the program 
so that the additional storage pained will allow any number 
of Q data points to be used. This will allow making the 
increments between data points smaller and extending the 
range, where possible, as follows: 


.001 £ N, < 1,000,000 
0001 <¢ < 1,000 


in order to make the approximation more accurate in the 
desired range of 


pO es N, = 000 2000 
.0O1 < gc < 100 
b) employ double precision arithmetic. 


c) permit using larger values of M, N, and P. 
3. Give more thought to the selection of appropriate functions 


of x/a, N. and @ which might improve the approximation without requiring 
2 


B 
larger values of M, N, and P. 
4. Once the values of peat and Y;s,% are evaluated, insure that 


they are properly placed in their respective matrices and the coefficients 


A determined correctly. 
mnp 
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5. Review the "Check Program" and insure that it correctly evalu- 


ates the approximate solution. All work should be done in double precision 


to insure the computer does not contribute to the error. 
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APPENDIX A 
HEAT TRANSFER CURVES 


Use of the curves is quite simple when the following procedure is 
applied: 
Ll. Select the appropriate curve. 
2. Determine the Biot number N, = ha/k. 
3. Determine the Fourier number © = Patt ae and logjg ©. 
4. Using log.) @ as abcissa, find the intersection with the 
appropriate curve corresponding to N, (interpolate Ne as necessary). 
5. Go to the left at this height on the graph. On the axis 
of ordinates read the value of the desired result. 


6. Perform any indicated arithmetic operation. 
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PLOT OF T(a)/F VS. LOG,9 OF FOURIER NUMBER FOR STEP TEMPERATURE CHANGE AT 
VARIOUS BLOT NUMBERS : 
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Fig. A-3 
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Fig. A-4 
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Fig. A-5_ 
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Fig. A-10 
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LOG) OF FOURIER NUMBER 
PLOT OF T(a)-Tavg/Et VS. LOG;9 OF FOURIER NUMBER FOR LINEAR TEMPERATURE 


CHANGE AT VARIOUS BIOT NUMBERS 


Fig. A-12 
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Fig. A-13 
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Fig. A-14 
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Section II 
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